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Abstract 

The coincidence problem is studied for the dark energy model of effective Yang-Mills condensate 
(YMC) in a flat expanding universe during the matter-dominated stage. The YMC energy Py{t) 
is taken to represent the dark energy, which is coupled either with the matter pm{t), or with both 
the matter and the radiation components Pr{t)- The effective YM Lagrangian is completely 
determined by quantum field theory up to 1-loop order with an energy scale k^/^ ~ lO^^ev 
as a model parameter, and for each coupling, there is an extra model parameter introduced. 
Beyond the non-coupling case, we have extensively studied four types of coupling models: 1) 
the YMC decaying into the matter at a rate V ~ 0.5i/o. where Hq is the Hubble constant; 
2) the matter decaying into the YMC at a rate V ~ 0.02ifo; 3) the YMC decaying into both 
the matter (at V ~ 0.5i/o) and the radiation (at V ~ 1.8 x lO^^iJo)- 4) both the matter (at 
r ~ O.Sifo) and the radiation (at V ~ 1.8 x IQ~^Hq) decaying into the YMC. In each of these 
four models, we have also explored various couplings. For all these models and for a variety of 
functional forms of the couplings F and V , it is found that the overall feature of the cosmic 
evolution for the YMC component is similar. Starting from the equality of radiation-matter 
Pmi = Pri, foi' generic initial conditions of pyi subdominant by a factor ranging over 8 orders 
of magnitude from 10^^° to 10~^, the models always have a scaling solution during the early 
stages, and the YMC always levels off at late time and becomes dominant, so that the universe 
transits from the matter-dominated into the dark energy-dominated stage at z ~ (0.3 ~ 0.5), 
and evolves to the present state with Vty ^ 0.7, Vtm — 0.3, and Vtr — 10^^. For the matter 
and radiation components, their evolutions depend on how they couple to the YMC. It is found 
that, if the YMC decays into a component with a rate being a constant or depending on py, as 
in Model 1 and Model 3, then in the early stages one has Pm{t) oc a{t)^^ and pr{t) oc a{t)~^ 
like the non-coupling case, but later around z ~ 0, this component also stops decreasing and 
levels off, asymptotically approaching to a constant value, just as py{t) does. In these cases, 
the equation of state (EoS) of the YMC Wy = Py/py crosses over —1 around 2; ~ 2 and takes 
the current value Wy ^ —1.1 at 2 = 0, consistent with the recent preliminary observations on 
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or if both the matter and radiation decay into the YMC in Model 4, then Pm{t) oc a{t)~^ and 
Pr{t) OC approximately for all the time, and Wy approaches to —1 but does not cross 

over —1. We have also demonstrated explicitly that, the coupled dynamics for {py{t), Pm{t)) 
in Model 1 , or for {py{t), Pm{t), Pr{t)) in Model 3, as i — > oo, is a stable attractor; in Model 
2 and Model 4 the dynamics of py{t) has a stable attractor as t — > oo. Therefore, under very 
generic circumstances, the existence of the scaling solution during the early stages and the 
subsequential exit from the scaling regime around z ~ (0.3 ~ 0.5) are inevitable. Thus the 
coincidence problem is naturally solved in the effective YMC dark energy models. 
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1. Introduction 



The observations on the cosmic microwave background radiation (CMB)[1] suggests a flat universe. This 
observation together with that from the Type la Supernova (SN la) [2] implies that the universe consists of 
some mysteries dark energy (f^A ~ 0.73), dark matter (fi^ ~ 0.23), and ordinary baryon matter (O5 ~ 0.04), 
and a radiation component Qr ~ 10~^. This is also supported by and the large scale structure of the universe 
[3] . The dark energy as the dominant cosmic energy component drives the current accelerating expansion of 
the universe. The simplest model for the dark energy is the cosmological constant A, which corresponds to a 
homogeneous and time-independent energy density pA ~ 5.8/io x 10~^^ev^ with an EoS of wa = Pa/ PA = — 1- 
Although this model can account for the observation so far, it has a coincidence problem. The observations 
show that the present value of the energy density for the matter component Pm = Pd + Pb is about one third 
of PA, but it varies with time as Pm{t) a{t)~^. So, for example, at an earlier time of radiation-matter 
equality with redshift z ~ 3454 [4], pA should be a very fine tuned value ~ 6.3 x 10~^^pjn{t). Otherwise, 
a slightly variant initial value of pA would lead to a value of the ratio Pa/ Pm drastically different from the 
observed one. This is called the coincidence problem. 

One class of models aiming at solving the coincidence problem is based upon the dynamics of some 
scalar field (/>, such as quintessence [5], K-essence [6], tachyon [7], phantom [8], and quintom [9], etc. These 
models can give rise to certain desired features of evolutional dynamics, such as scaling solutions [10] and 
tracking behavior. As a common point, these scalar models need to make use of some special forms of the 
potential V{(f)) with certain chosen parameters. For example, in the quintessence model, one needs to choose 
V{(f)) oc {M/(j))" with a and M being a positive number, or oc e^*/"'?'' [11]. In phantom model one may take 
y((/)) = Vo[cosh(Q(/)/mp/)]^^ [12]. Moreover, phantom models typically introduce a negative kinetic energy 
term —(fP'/2. Sonic of these models are expected to be low energy effective field theory, coming from some 
fundamental field theory, others are simply introduced by hand. 

As far as cosmological observations are concerned, a cosmological model has to give the current status: 
Oa ~ 0.7 and Q„( ^ 0.3. Besides, it is safer for a model not to contradict the conventional scenario in the Big 
Bang model from the energy scale ~ IMev down to the present. Therefore, one idea to solve the coincidence 
problem is that during the early stages of the expansion the dark energy density need not to be a constant, 
but varies with time. Even its EoS w need not to be close to —1 in early stages. Only at some rather 
recent moment has it become dominant and acquired an EoS u; ~ — 1. In order to allow the conventional 
cosmological processes, such as the nucleosynthesis and the recombination, etc., to have been occurred in the 
past, the dark energy component should be subdominant to the matter component early stages. 

In the approach of the scalar field models, certain coupling has been introduced between the scalar field 
dark energy and the matter [13] [14]. With some particular choice of the model parameters, there can exist 
a scaling solution of dynamics, in which the dark energy density is proportional to that of matter during the 
early stages of evolution. However, to achieve the tracking solution, some scalar models need to have a very 
large coupling, so that the universe would enter the acceleration stage soon after the matter era. This would 
result in a picture of structure formation, totally different from that required by the observations. To remedy 
this defect some models [15] [16] introduce certain particular form of couplings, but still the entrance to the 
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out that the k-essence models always have, at some stage, the difficulty of superluminal propagation, leading 
to violation of causality, and the EoS and the sound speed of k-essence could be greater than > 1 [17]. Besides 
the unconventional negative kinetic energy, after the EoS w cross —1, a class of scalar models may suffer 
from severe quantum instabilities and the Big Rip singularity, i.e. either the energy density p, or the pressure 
p, may grow to infinity within a finite time. More recently, an overall estimate of scalar models has been 
given on the issue of coincidence problem. By examining the most general form of scalar Lagrangian with a 
generic coupling between scalar dark energy and dark matter, it has been shown that a vast class of scalar 
models, including all the models in current literature, have the difficulty to implement both a scaling solution 
without singularity and a sequence of expansion epochs that required by standard cosmology, such as the 
radiation-dominated, matter-dominated, and dark energy-dominated epochs [18]. Therefore, the coincidence 
problem still remains after over a decade of extensive studies [19]. 

The introduction of the quantum effective YMC into cosmology [20] has been motivated by the fact that 
the SU (3) YMC has given a phcnomcnological description of the vacuum within hadrons confining quarks, 
and yet at the same time all the important properties of a proper quantum field are kept, such as the Lorentz 
invariance, the gauge symmetry, and the correct trace anomaly [21]. Quarks inside a hadron would experience 
the existence of the Bag constant, S, which is equivalent to an energy density p = B and a pressure p = —B. 
So quarks would feel an energy-momentum tensor of the vacuum as Tj^^, — B dicLg(\, — 1, — 1, — 1). This non- 
trivial vacuum has been formed mainly by the contributions from the quantum effective YMC, and from the 
possible interactions with quarks. Our thinking has been that, like the vacuum of QCD inside a hadron, 
what if the vacuum of the universe as a whole is also filled with some kind of YMC. Gauge fields play a 
very important role in, and are the indispensable cornerstone to, particle physics. All known fundamental 
interactions between particles are mediated through gauge bosons. Generally speaking, as a gauge field, 
the YMC under consideration may have interactions with other species of particles in the universe. In our 
previous studies of [23, 24], the possible interaction of the YMC with other cosmic components have not been 
examined. However, unlike those well known interactions in QED, QCD, and the electro-weak unification, 
here at the moment we do not yet have a model for the details of the microscopic interactions between 
the YMC and other particle. Therefore, in this paper on the dark energy model, we will adopt a simple 
description of the possible interactions between the YMC and other cosmic particles. That is, we introduce 
coupling terms in the continuity equations of cosmic energy densities, such as the YMC, the matter and the 
radiation, study the cosmic evolution of the universe from the matter dominated era up to the present. As 
shall be seen, in our model the current status of the universe turns out to be a natural result of evolutional 
dynamics driven by the effective YMC as the dark energy, plus the matter and the radiation that are coupled 
to the former. What is important is that this has been achieved with a choice of the initial value of the 
fractional energy density of the YMC ranging from 10~^° to 10~^. As a novel feature in contrast to the 
non-coupling models, the coupling YMC dark energy models can give rise to an EoS of dark energy Wy 
crossing —1, say 'u;^ ~ — 1.1 at present. The coincidence problem can be solved. 

In section 2, as a basis for the setup, an introduction is given to the effective Yang-Mills condensate 
theory, and the dynamic equations for the three components in the Robertson- Walker spacetime are derived. 
Section 3 is about the simplest case of non-coupling. Section 4 studies the dynamic cosmic evolution with 



the effective YMC decaying into tfie matter component. Tfiere are two dynamic equations for the YMC 
and the matter, respectively. Section 5 studies the model in which the matter decays into the YMC. The 
matter component has different dynamic evolutions especially at late stages for these two coupling models. 
Nevertheless, in both these two models, the YMC has a scaling solution and a tracking behavior as a natural 
outcome of the dynamic evolution. Section 6 extends to the general case with the YMC coupling to both 
the matter and the radiation. Now one has one more coupled dynamic equation and an extra coupling for 
the radiation component. Two cases are studied, the YMC decaying into the matter and radiation, and the 
matter and radiation decaying into the YMC. The dynamics is examined parallel to sections 4 and 5. In 
each of these four models, various functional forms of coupling have been explored. The major part of the 
study is in sections 4, 5, and 6, which contain the calculations and the results. Section 7 gives an analysis 
on asymptotic behavior of the dynamic system at t ^ cx). It is found that there exists a unique attractor in 
the asymptotic region, which is stable against perturbations. Section 8 contains a summary and discussions. 
Throughout this paper we will work with unit, in which c = ?i. = fcs = 1. 

2. YM condensate as dark energy 

In the effective YMC dark energy model, the effective YM field Lagrangian is given by [20] [21]: 
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where k is the renormalization scale of dimension of squared mass, F = -\F^^F°■^''' = E"^ - B'^ plays the 
role of the order parameter of the YMC. In this paper, for simplicity, we only discuss the pure 'electric' case, 
F = E"^. The Callan-Symanzik coefficient b = (UN - 2Nf)/247r'^ for SU{N) with Nf being the number 
of quark flavors. For the gauge group SU{2) considered in this paper, one has 6 = 2- ll/247r2 when the 
fermion's contribution is neglected, and 6 = 2 - 5/247r2 when the number of quark flavors is taken to be 
Nf = 6. For the case of SU{?)) the effective Lagrangian in Eq.(l) leads to a phenomenological description of 
the asymptotic freedom for the quarks inside hadrons [21]. It should be noticed that the SU{2) YM field is 
introduced here as a model for the cosmic dark energy, it may not be directly identified as the QCD gluon 
fields, nor the weak-electromagnetic unification gauge fields, such as Z° and W^. As will be seen later, 
the YMC has an energy scale characterized by the parameter k^/^ ~ 10~^ ev, much smaller than that of 
QCD and of the weak-electromagnetic unification. An explanation can be given for the form in Eq.(l) as an 
effective Lagrangian up to 1-loop quantum correction [21]. A classical SU{N) YM field Lagrangian is 
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where go is the bare coupling constant. As is known, when the 1-loop quantum corrections are included, the 
bare coupling constant go will be replaced by the running coupling g as the following [21] [22] 
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where k is the momentum transfer and ko is the energy scale. To build up an effective theory [21], one may 
just replace the momentum transfer /e^ by the field strength F in the following manner: 
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yielding Eq.(l). We like to point out that the renormalization scale k is the only parameter of this effective 
YM model, and its value should be determined by comparing the observations. In contrast to the scalar-field 
dark energy models, the YMC Lagrangian is completely fixed by quantum corrections up to order of 1-loops, 
and there is no room for adjusting its functional form. This is an attractive feature of the effective YMC 
dark energy model. 

Prom Eq.(l) we can derive the energy density and the pressure of the condensate [20, 23] in the flat R-W 
spacetime: 

Py = \^E^ + \bE\ (2) 
Py = \eE^-\bE\ (3) 
where e is called the dielectric constant of the YMC, given by [21] 
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MnK|. (4) 



The EoS of YMC is given by 



where 



Py 3y + 3 



is a dimensionless quantity, in terms of which the energy density and pressure of the YMC will be given by 

Py = lbK\y + l)ey, (7) 

Py = hK\^y-l)ey. (8) 

One sees that, to ensure that the energy density be positive in any physically viable model, the allowance 
for the quantity y should be y > —1, i.e. F > K^/e ~ 0.368k^. Before setting up a cosmological model, 
the EoS Wy itself as a function of F is interesting. From Eqs.(2) and (3) one sees that the YMC exhibits 
an EoS of radiation with py = and tt; = 1/3 for a large dielectric e ^ b (i.e. F ^ k^). On the other 
hand, for e = (i.e. F = kP), which is called the critical point, the YMC has an EoS of the cosmological 
constant with py = —py and w = —1. The latter case occurs when the YMC energy density takes on the 
value of the critical energy density py = ^bK^ [20]. It is this interesting property of the EoS of YMC, going 
from w = 1/3 at higher energies (F S> k^) to if; = — 1 at low energies {F = k^), that makes it possible 
for the scaling solution [10] for the dark energy component to exist in our model. More interestingly, this 
transition is smooth since w is smooth function of y in the range (—1, oo). Now we ask the question: Can 
Wy cross over —1? By looking at Eq.(5) for Wy, we see that Wy only depends on the value of the condensate 
strength F. In principle, u;^ < — 1 can be achieved as soon as F < k^. Moreover, in regards to the behavior 
of tWj, as a function of F, this crossing is also smooth. However, as shall be shown explicitly later, when 
the YMC is put into a cosmological model as the dark energy component, together with the radiation and 
matter components, to drive the expansion of the universe, the value of F can not be arbitrary, it comes out 
as a function of time t and has to be determined by the dynamic evolution. Specifically, when the YMC does 
not decav into the matter and radiation, w„ can onlv anoroaches to —1 asvmntoticallv, but will not cross 



over —1. On the other hand, when the YMC decays into the matter and/or radiation, Wy does cross over 
—1, and, depending on the strength of the couphng, Wy will settle down to an asymptotic value ~ —1.17. As 
a merit, in this lower region oi Wy < —1, all the physical quantities Py, Py, and Wy behave smoothly, there is 
no finite-time singularities that are suffered by a class of scalar models. 

Now we put the YMC in to the cosmic setting, which is assumed to be a spatially fiat {k = 0) Robertson- 
Walker spacetime 

= _ a'^{t)8.jdx'dxK (9) 

As it stands, the present universe is filled with three kinds of major energy components, the dark energy, 
the matter, including both baryons and dark matter, and the radiation. In our model, the dark energy 
component is represented by the YMC, and the matter component is simply described by a non-relativistic 
dust with negligible pressure, and the radiation component consists of CMB and possibly other particles, 
such as neutrinos, if they are massless. Since the universe is assumed to be flat, the sum of the fraction 
densities is O = -|- $7^ + = 1, where the fractional energy densities are 0^ = Py/ p, = Pm/ P^ and 
Or = Pr/ P- The overall expansion of the universe is determined by the Priedmann equations 
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in which all these three components of energy contribute to the source on the right-hand side of the equations. 
The dynamical evolutions of the three components are determined by their equations of of motion, which 
can be written as equations of conservation of energy [20, 23]: 

Py + ^'^{Py+Py) = -Qm-Qr, (12) 
Pm + = Qm, (13) 

a 

Pr + ^'^{Pr +Vr) = Qr, (14) 

where Qm represents the energy exchange between the YMC and the matter, and Qy. between the YMC 
and the radiation, respectively. In the natural unit, both quantities have the dimension of [energy]^. The 
couplings Qm and Qr are phenomenological, and their specific forms of will be addressed later. The sum 
of Eqs. (12), (13), and (14) guarantees that the total energy is still conserved. As is known, Eq.(ll) is not 
independent and can be derived from Eqs. (10), (12), (13), and (14). It is noted that once the couplings Qm 
and Qr are introduced as above, they will bring two new parameters in our model. When Q^ > 0, the YMC 
transfers energy into the matter, and this could be implemented, for instance, by the processes with the 
YMC decaying into pairs of matter particles. On the other hand, when Qm. < 0, the matter transfers energy 
into the YM condensate. Similarly, when Q^ > 0, the YMC transfers energy into the radiation. Therefore, 
in the most general case of coupling, there will be three model parameters: Qm, Qr, and k. 

In the following computations, it is simpler to employ the following functions rescaled by the critical 
energy density ^bK^ of YMC, 



Qm rT7\ 

= w ' * 

Here the dimensionless functions x and r are simply the rescaled energy density of the matter and radiation, 
respectively, and the rescaled exchange rates Qm and qr have unit of [time]~^. Then, in terms of x, y, and r, 
the dynamical evolutions given in Eqs. (10), (12)-(14) can be recast into: 

dy , 4y q,n + qr 



dN 2 + y Hh{2 + y)ey' 



(19) 



{^f = H^h^ (22) 



where the variable A'' = Ina(i), the function h = -\/(l + y)ey + x + r, and the constant H = ^/AirGbK^/S. 
Note that H is not exactly the present Hubble constant Hq. Prom Eq.(22) one can see that it is the quantity 
Hh that determines the actual expansion rate of the universe, and the present value of Hh is identified as 
the Hubble constant Hq. However, as our calculations will show later, the value of /i ~ 1.07 at present, so 
approximately H Hq. It should be emphasized that H is not an independent parameter, it is fixed by 
the model parameter k. Once the YMC is put into the cosmological context, one can estimate the order of 
magnitude of k as follows. The critical density pc = 8.099/io x lO^^^ev^ with the Hubble parameter ho 2± 0.72, 
and the current value of the dark energy density should be py = ClyPc c± 0.7pc- As calculations will show, the 
present value of the factor (1 + y)e^ ~ 0.8 in Eq.(7), and py 2± 0.8 x ^6ac^. So one has 

~ 5 X IQ-^ho ev. (23) 

This energy scale is much smaller than those typical energy scales occurring in the standard model of particle 
physics, such as ~ 10^ Mev for QCD, and ~ 10^ Gev for the weak-electromagnetic unification. Therefore, for 
lack of an explanation of the origin of this energy scale k within the standard model, we may have to regard 
the YM condensate as a new physics beyond the standard model of particle physics. In this sense, like other 
dark energy models, the fine-tuning problem, i.e. why k has just this small value, also exists in our model. 

The set of equations (19) - (22) hold for a generic stage of cosmic expansion driven by the combination 
of the radiation, the matter, and the dark energy. In this paper we focus on the matter-dominated era and 
the subsequent accelerating era. In particular, we like to see how the cosmic expansion evolves and transits 
from the matter-dominated to the accelerating era. In the remaining of the paper, we always take the initial 
condition to be at the time U of the equality of radiation-matter (with a redshift z = 3454l3g2 [4]) 



(24) 



with the subscript i denoting the initial value. Of course, once the initial value flyi for the YMC is given, 
one has immediately 
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In order to keep the main features of the conventional scenario of cosmic expansion, it is also assumed that 
initially the matter and the radiation are dominant, = ^ri — 1/2, and the YMC is subdominant 

< 1/2. (26) 

In Ref.[24] we considered the constraints on the YMC energy density at an earlier stage. There the initial 
condition was taken at a redshift z ~ lO^'^, corresponding to an energy ~ 1 Mev, during the radiation 
stage, when the Big Bang nucleosynthesis processes took place. The upper bound has been found to be 
Uy < 0.26 at z ~ lO^'^. Afterwards up to z ~ 3000, both pr and py evolved approximately in a similar 
way oc a{t)~^. Thus, at the equality of radiation-matter with z 2± 3500, this consideration will give an upper 
bound Q,yi ~ 0.1. Within this restriction, the initial value flyi may still be allowed to vary in a very broad 
range. In the following we take a safe value 

< 10"^ (27) 

as the upper bound both for illustration purpose. This choice has been made in concordance with the thinking 
that the dark energy component has been existing in the universe from the equality of radiation- matter, but 
its initial relative contribution is bounded by a few percent of the total, so that during most of the history of 
the universe the cosmic evolution will be the same as in the standard Big Bang model. Only quite recently 
has the dark energy component become dominant and modified the cosmic evolution considerably. 

3. Non-coupling Case 

The simplest case is the non-coupling with Q„i = and Qr = in Eqs.(12), (13), and (14). Then there is 
only one model parameter k that has already been fixed in Eq.(23). Each component evolves independently 
in the expanding RW spacetime. To solve the dynamic equations, we take the initial YMC at the time ti to 
be in the broad range 

Qyi = (10-^0, 10-2), (28) 

consistent with the restriction (27), that is, the initial value of pyi ranges over eight orders of magnitude. In 
terms of y, x, and r, the initial condition above can be written as 

Vi = (1, 16.12), (29) 

Xi = ri = 1.7 X 10^°. (30) 

The equations (19) (20) (21) with Qm = Qr = ^ are solved easily for Py{f), pm{t), and Pr{t), which are shown 
as a function of redshift z in Fig.l. Both pjn{t) oc a~^{t) and pr{t) oc a~'^{t) decrease monotonically at their 
fixed slope, respectively, and do not level oS as t ^ oo. More interesting is the evolution of YMC energy 
density. In the early stage py{t) is subdominant to Pm{t) and Pr{t), and decreases at a slope between those 
of pr{t) and Pm{t), tracking the matter. Later, Py{t) gradually levels off at and approaches to a constant. 



At z ~ 0.35, Py{t) starts to dominate over Pmit) and the accelerating expansion takes over. This exit from 
the subdominant region (scahng) to the dominant region is naturally realized. It is important to notice that, 
as long as the initial value is in the broad range of Eq.(28), Py{t) always has the same asymptotic value 
as t — oo. By the way, the first order differential equations (20) and (21) for x{t) and r{t) have no fixed 
points since dx/dN ^ and dr/dN ^ during the course of evolution, but the differential equation (19) for 
y{t) has a fixed point = as solution of dy/dN = as t — oo. Fig.2 gives the corresponding evolution 
of the fractional energy densities. Starting with the initial value Vtri 1/2, the radiation component VLj. 
has a simple evolution of monotonic decrease. In contrast, the matter component $1^, starting with c^i 1/2, 
increases quickly and approaches to ~ 1 around a redshift (1 + z) ~ 174. At (1 + z) ~ 2.7, $7^ drops down 
and is dominated by Vly at z ^ 0.35. The YMC component Q.y starts with the very small initial value and 
increases slowly and monotonically. Around (1 + z) ~ 2.7, Vly has a quick increase, and around z ~ 0.35 it 
dominates over Q^n- Observe in Fig.2 that the two curves of Qy for the two different initial values lO""^'^ and 
10~^, respectively, are almost overlapped into one curve. This pattern of degeneracy for Q,y demonstrates 
vividly the fact that the cosmic evolution and the current status are insensitive to the initial condition of 
the YMC. As the result of evolution, at present (z = 0), one has Q.y ~ 0.7, Q,m ^ 0.3, and 9,r - 10"^ Fig.3 
shows the evolution of Wy and that of the effective EoS Weff defined by 

Both approach to —1 as t ^ oo, but, they can not across —1. So there is no super-accelerating stage in 
the non-coupling model. Looking back Eq.(5) and (6), it is clear that in the non-coupling case the YM field 
strength will always stay above the critical value: F > k^. The asymptotic region at t — oo corresponds to 
F = k!^. This is a state with the dielectric constant e = 0. Thus, for the non-coupling case, no matter what 
kind of initial condition is given, the YMC always settles down to the state of e = as a result of dynamic 
evolution. As for the current status Cly ~ 0.7 and fi^ ~ 0.3, it has been achieved for the whole range of 
initial values in Eq.(29) at the fixed model parameter k in Eq.(23). Therefore, the coincidence problem is 
solved in this model, but the fine-tuning problem still exists, i.e., wc do not have an answer to the question 
why K should have such a value as in Eq.(23). The non-coupling case with Qm = and Qr = has also been 
studied in [23, 24] with the initial condition being taken at a redshift z ~ 10^*^ in the radiation dominated 
stage. There the evolution behavior found for the matter dominated era is similar to what is obtained here. 

4. YMC Decaying into Matter 

Consider the case that the YMC couples to the matter component only. In terms of the cosmic energy 
densities today, the radiation fraction is roughly $7^ ~ 10~^, a relatively very small contribution, much lower 
than the other two components. Therefore, in this section we temporarily neglect its coupling with the YMC 
by setting Qj. = 0. There is only one free parameter Q^n 

since k has been fixed. Then Eqs.(19), (20), and 

(21) reduce to 

dy 'iy _ Qm .o2^ 

dN 2 + y Hh{2 + y)ey ' ^ ^ 

^ + ^- = ^^ (33) 



- + 4. = 0. (34) 

The dynamic evolution of the radiation is independent of the other two, and pr{t) oc a^^(t), as is seen in 
Eq.(34). To proceed further, one needs to know the coupling qm to solve Eqs. (32) and (33). As mentioned 
earlier, in general there are two kinds of models depending on whether Q„i > or Q„i < 0. In in section, we 
examine the model Qm > with the YMC decaying constantly into the matter. We call this the Model 1. 
It can be generically expressed as 

Qm = ^Py, (35) 

where F > is of dimension [time]~^ and measures the decay rate of the YMC energy density as well as 
the production rate of the matter energy density. Here T is simply taken as a model parameter describing 
phenomenologically the interactions between the YMC and the matter. In the following we discuss several 
cases for the parameter T. Substituting Eq.(35) into Eqs. (32) and (33) yields 

dx r (1 + y)ey 



dN H h 
dy T l + y Ay 



- 3x, (36) 

(37) 



dN H{2 + y)h 2 + y 
1. Consider the simple case of a constant rate with 

r/H = 0.5. (38) 

We have taken this magnitude of the rate F, so that the present status of the universe will be Qy ~ 0.7 and 
fljn ~ 0-3 as the outcome from our computation. Interestingly, in order to achieve this status, as a model 
parameter, the decay rate of the YMC needs to be of the same order of magnitude as the expansion rate of 
the universe, i.e. F H. Now at the time U with z ~ 3454 the initial YMC energy density is taken to be in 
the range 

Qyi = (10-^°, 3 X 10"^), (39) 

similar to Eq.(28), which corresponds to 

Vi = (1, 15). (40) 
The initial values for the matter and the radiation components are given by 

Xi = ri = 1.0x10^^, (41) 

which are a little bit smaller than that in Eq.(30). This is because in the case here the matter are being 
generated out of the decaying YMC during the course of evolution. Consequently, smaller initial values Xj 
and Vi are needed to arrive at the current status. 

Given the different initial values of Clyi in Eq.(39), the corresponding initial values of Jl^j = Uj-i = 
(1 — Q.yi)/2 also vary by a small amount. But for the different values of yi in Eq.(40) we have taken the 
same set of values in Eq.(41) since the initial total energy density pyi + pmi + pn itself at z ~ 3454 has some 
errors. The results are given as functions of the redshift z in Fig.4, Fig.5, and Fig.6 . The evolution of Py{t) 
is also similar to the non-coupling case. During the early stage Py{t) is lower than, and keeps track of pmit)- 
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and the accelerating stage begins. In fact in Fig.4 for three different yi there are three curves of the matter 
Pm{t), respectively. However, these three curves are too close to each other so that they are overlapped. The 
similar is for Pr{t). Note that, in the presence of coupling Q^^, the evolution of Pm{t) is different from the 
non-coupling case, in that at the late stage around z ~ it also levels off just like Py{t). In this sense, the 
evolution of the matter component is sort of bound to the YMC. In fact, as t ^ oo, the set of Eqs.(36) and 
(37) have the asymptotic behavior 

and both py{t) and pm{t) have asymptotic values. This will be addressed later in section 6. 

As a novel feature of the coupling model, in contrast to the non-coupling case, now the EoS of the YMC 
as a function of time t, Wy{t) crosses over —1 around z ^ 2, takes a value Wy ~ —1.1 at present ^ = 0, and 
approaches u;^ ~ — 1.17 asymptotically, as shown in Fig.6. This occurrence of crossing over —1 in this model 
can be understood, since the coupling makes the YMC to loss energy into the matter, consequently, the YM 
field strength F will drop down below the critical value k^, leading to e = 6y < and Wy < —\ as \n Eq.(5). 
This can also be arrived by looking at the asymptotic region determined by the equation ^ = 0, which by 
Eq.(37) is just 

ii±i^+% = 0. (43) 

Recall that for the non-coupling F = the asymptotic value is y/ = 0, yielding Wy = —1 hy Eq.(5). Once 
F > 0, Eq.(43) yields an asymptotic value < as the solution, hence Wy < —1. Thus, when transferring 
energy to the matter, the YMC will eventually settle down in the state of Wy < —1, which is equivalent to a 
negative dielectric e < 0. Recently, there are some observational indications that the current value of EoS of 
dark energy w is less than —1, for instance, w = —1.023 it 0.090{stat) it 0.054(sys) from the 71 high redshift 
supernovae discovered during the first year SNLS [25], and w = — 1.2lto}2 from the blind analysis of 21 
high redshift supernovae by CMAGIC technique [26]. Of course, this is still to be observationally examined 
with higher confidence level in future. However, the crossing over —1 would be difficult for scalar models, 
except for quintom models at a price of introducing two scalar fields and an artificially designed potential 
[27]. As we just have demonstrated, in the YMC dark energy model with coupling > 0, this crossing 
is realized naturally. In general, the asymptotic value of u;^ at f — oo and the current value Wy Sit z = 
as well, are determined by the asymptotic value of y through Eq.(5), and the latter is obtained from the 
combination of Eqs.(37) and (42), and thus depends on the ratio F/H. Thus the asymptotic value of Wy is 
determined by the ratio F/H of the two parameters of our model. For instance, the observed Eos of dark 
energy w = —1.023 from SNLS [25] can be obtained, in our model, by taking a little smaller decay rate 
F/H = 0.13 and sHghtly higher initial densities = = 1.5 x 10-'^'^, and the value w = —1.21 from the 
analysis by CMAGIC [26] can be obtained by taking a bit larger decay rate F/H = 0.81 and a slightly lower 
densities Xi = Vi = 0.7 x 10^*^, respectively. Fig.6 also shows that the effective Weff can not cross —1 yet, 
and its asymptotic value is ~ —0.96. Interestingly, this model predicts that in the upcoming future the dark 
energy density py will remain a constant slightly larger than today, and that the matter energy density p^Yi 
will be a constant slightly lower than today. Eventually the universe will settle down to a steady state with 
9.y ~ 0.85, 0.15, and - 0. 



Therefore, this model yields a picture of evolutional cosmos, the early part of which can account for the 
the past history of the expanding universe, i.e., that of the standard Big-Bang model, and the late part of 
which, i.e. the future of the universe, is similar to the Steady State Model [28] [29]. As a plus for the YMC 
model, there are no Big Rip singularities in finite time, since all the quantities py, py, and Wy are smooth 
function of the time t. 

Although in the above we have only presented the results for the initial values within the range Q,yi > 10~^^ 
given in Eq.(28), as a matter of fact, we have worked out for even lower values ^yi < 10"-'^'^. For instance, 
we have calculated the case of ^yi ~ 2 x 10~^^, i.e., yi = —0.9. the initial energy density pyi are even lower 
than the current values py ~ 0.7pc- The details are given by the curves denoted by yi = —0.9 in Fig.4, Fig.5, 
and Fig. 6. The evolution is such that py{t) starts initially from the given very low value, increases (instead 
of decreasing) very quickly with time t , and approaches its corresponding asymptotic value ~ 0.7 Pc, as is 
seen in Fig.4. The evolutions for Pm{i) and pr{t) arc similar to the cases in Eq.(39). Thus, except the initial 
increasing of Py{t), the current status of the universe is the same as those in the range of Eq.(39). 

2. Now consider a case of the YMC decaying into fermion pairs. As is known in QED, a constant 
electric field is unstable against decay into pair of particles [30]. Analogously, a constant 'electric' SU{3) 
Yang-Mills field configuration of QCD is unstable against decay into quark pairs and gluon pairs [31] [32] 
[33]. Calculations have shown that for a gauge group SU{N), due to the decay into fermions pairs and gauge 
boson pairs, the average local color electric field decreases at a rate 

^ = -cE^^ (44) 

where c is a dimensionless constant [32] [33] [20] . Notice that this rate is derived from the energy- momentum 
conversation of pair creation. From the point of view of particle physics, the SU{2) YM field in our case 
should be allowed to have some fundamental interactions with other microscopic particles. At the moment 
we do not intend to proceed further to build up the detail of its interaction. Instead, we simply assume that, 
similar to the SU{3) gauge field in QCD, the YM condensate, here a constant in space, also has the property 
of being unstable against decay into other particles, and the decay rate of the field strength has the same 
form as Eq. (44) . Consequently, this in turn will give rise to the decay rate of the YMC energy density 

Making use of the chain rule relation ^ = ^ ^ and Eq. (44) , one obtains the following expression for the 
decay rate [20] 

r = 2cK^^^e4. (46) 

depending on the coefficient constant c, which is treated as the model parameter in place of F. We choose 
the value of c = O-UbH/K^/"^ , so that 

^ = 0.25^ef. (47) 

H l+y ^ ' 

As our computations show, the value of variable y at present stage is very small y ~ 0, so Eq.(47) yields 
F ~ O.Sif, quite close to that in Eq.(38). The initial condition is taken to be the same as Eqs.(40) and (41). 
Figs 7, 8, and 9 show the results, which are very similar to that in the previous case of the constant rate. 
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estimated as ~ ^Hp^^ — g cm~^ , a very small rate, equivalent to generation of ~ 0.3 protons in 
a cubic kilometer per year. This value is approximately equal to that in the Steady State Model [28] [29]. 
Thus in our model the particle pairs are continuously generated, at a very low rate, out of the vacuum filled 
with the YMC. 

As is known, in order to continuously generate the cosmic matter, the Steady State Model has to introduce 
some C-field with negative energy [29], which is problematic in a physical theory. Here in our model, it is 
the effective quantum YMC plays the role of a matter-generator, there is no negative energy to occur in the 
proper range. The YMC has a positive energy and a negative pressure. 

Prom Eq. (38) and (47) it seems that the overall behavior of the dynamic evolution is not sensitive to the 
particular form of the coupling F, as long as its magnitude is F 0.5H. This has been confirmed in our 
examinations. For example, we have also investigated the case with the decay rate of the form 

T/H = 0.5e-y, (48) 

and the results are very similar to that in the previous case. Therefore, in Model 1 of the YM-matter coupling 
with the coupling = Fpy, as long as F is constant or depends on py, the overall features of the dynamic 
evolution are similar. We have also studied the cases with F depending on the matter p^. It is found that, 
if the decay rate F depends on the matter pm, for instance 

T/H = bx, (49) 

where b is some constant and x is defined in Eq.(15), then for b < 10"^ the evolution will be similar to the 
non-coupling case. When the constant b > 10~^, the decay rate is too fast, and at start Py{t) drops down 
quickly, later it increases to its asymptotic value from below. To keep the paper short, we do not demonstrate 
these detailed graphs of the cases of Eqs.(48) and (49). 

5. Matter Decaying into YMC 

In section we study the situation in which the matter decaying constantly into the YMC with < 0, 
just opposite to the Model 1. We call this the Model 2. It can be generically expressed as 

Qm = -Tpm, (50) 

i.e. the matter transports energy into the YM condensate. Then Eqs. (32) and (33) reduce to 

dx T X 

— + 3x = ---, 

dy Ay r X 
dN'^ 2 + y^ H{2 + y)eyh' 

1. Consider the simple case of a constant decay rate 

T/H = 0.02. (53) 

Again, here the value 0.02 has been taken, so that the resulting energy densities from our computation will 
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(51) 
(52) 



needs to be almost two order of magnitude smaller than the expansion rate. The initial values of Xi are taken 
to be 

Xi = riC^ 1.8 X 10^°. (54) 

Since the matter is decaying and constantly being converted into the YMC component, a slightly larger 
initial value of the matter has been taken than that in Eq.(41). The initial value of the YMC is taken to 
be Ui = (1,15), the same as Eq.(40). The results are given in Figs. 10, 11, 12. Now Wy does not cross over 
— 1. It is interesting to find out that the evolution of the YM condensate behaves differently for two different 
ranges of j/j. For the higher range 

Vi = (5, 15), (55) 

corresponding to Clyi = (5 x 10~^, 3 x 10~^), all the quantities have similar evolution as the non-coupling 
case in Figs. 1, 2, 3. For the lower range 

Vi = (1, 5), (56) 

however, the YMC has an instantly sudden increase during the initial stage, and quickly catches up the 
evolution pattern of the j/j = 5 case. This is in contrast to the smooth behavior on the higher range. Thus, 
in order to have a rather smooth evolution within the Model 2, the initial value yi should be given by the 
higher range in Eq.(55). Moreover, we have also found that the overall behavior of the dynamic evolution is 
not sensitive to the particular form of the coupling Q^. For instance, we have checked a case of 

r/H = 0.02e"^, (57) 

and the resulting evolutions are similar to the case of Eq.(53). Thus, the coincidence problem can also be 
solved in Model 2. 

So far in Model 1 and Model 2, in regards to the coupling between the YMC and the the matter, we have 
not explicitly distinguished the baryons and the dark matter, and have assumed, for simplicity, the same 
coupling Qm for both the baryons and the dark matter. We can roughly estimate the current value of the cross 
section corresponding to the collisions involving the baryons as in Eqs.(38), (47), (48), (49) with b < 10~^, 
(53), and (57). For instance, take the baryon decay rate T hxH as in Eq.(49). Then, by definition, the 
rate is F ~ vna, where v is the baryon velocity, the baryon number density is n = pb/mh ~ 0.04/9c/m6, and 
a is the crossing section for the collisions between the baryons and the YM gauge bosons for this type of 
interaction . Then we can get an estimate: 

a ~ 25bmbH/vpc. (58) 

Taking b ~ 10"^, v ~ lO^km/s, the baryon mass mb ~ 0.94 Gev, the current Hubble constant for H, one 
has cr ~ 6 X 10~^^ cm^. This is an order lower than the Thomson's cross section ar — 6.7 x 10~^^ cm^ for 
in QED. Similarly, letting F ~ 0.5H as in Eqs.(38) and (47) for the YMC decaying into baryons, we would 
get (7 ~ 20 X 10~^^ cm^ analogously, slightly greater than ar- Therefore, given this magnitude for the cross 
section a in both cases, we would, in principle, be able to observe this kind of interactions occurring, either 
with the baryon being decaying into the YM boson pairs or the baryon pairs jumping out of the vacuum. 
However, as said earlier, the rate F for this kind of events is too low, giving Qm ~ 0-3 proton generated in 



one cubic kilometer per year. For the Galaxy of a volume ~ lO^(kpc)^, this rate is roughly equivalent to 
an amount of mass ~ lO^^M© generated per year, a small production rate. The chance may be small for 
directly detecting the event. Even if future experiments rule out or restrict the coupling with the baryons, 
one has to drop it or reduce its magnitude as a model parameter. Nevertheless, the coupling with the dark 
matter probably still remains. This is because the dark matter is usually assumed not to have interactions 
with ordinary particles, such as baryons, photons, etc. So it is difficult to directly detect productions of dark 
particle pairs and decays of dark particles. 

6. Coupling with Both Matter and Radiation 

It is quite natural to allow the YMC to couple with both the matter and the radiation simultaneously. 
Now we study Model 3 that the YMC decays into the matter and the radiation as well: 



Then Eqs. (19) - (21) reduce to 



Qm = ^Py > 0, Qr = T' Py > 0. (59) 
dy r + r' 1 + y 4y 
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dN ~ H h 

Consider the case of the constant decay rates 



- 3x, (61) 

- 4r. (62) 



r/H = 0.5, T'/H = 1.8 X 10"''. (63) 

Note that T' is lower than F by three orders of magnitude. These values of coupling are taken so that the 
current values are ^ly ~ 0.7, — 0.3, ~ 8.6 x 10~^ (including massless neutrinos). The initial condition 
is the same as in Eqs. (40) and (41). The results are given in Figs. 13, 14, 15. As before, the particular form 
of the couplings is not important for the overall behavior of evolution. For instance, we have also examined 
the case 

r/i/ = 0.25^ei, F'/ii- = 0.9 X 10-4 ^ei, (64) 

1 + y 1 + y 

based on an analogous consideration to Eq.(47). The evolution is similar to the case of Eq.(63). In these two 
cases of Model 3, due to the couplings with the YMC, both the energy densities, pm and pr, level off around 
z ~ 0, and Wy crosses over —1 around z ~ 2.5. Thus all the three components of cosmic energy will remain 
constant in future, and the state of the universe will keep almost as it is today. This is similar to the Steady 
State universe [28] [29]. Hence, according to Model 3, the past history of the universe is consistent with the 
conventional standard Big Bang model, and from now on, the cosmic evolution tends to that of a steady 
state, that is, the universe will remain similar to it is today, but with Qy ~ 0.85, Qm 0.15, and Qr W~^. 

We have also studied the case with both the matter and radiation decaying into the YMC. We call this 
Model 4. The couplings are such that 

Qm = -rpm < 0, Qr = -T' Pr < 0. (65) 



We take 

V/H = 0.02, T'/H = 1.8 X 10"^. (66) 

The initial condition for the YMC is the same as Eqs.(55) and (56) and the initial densities for the matter 
and radiation are the same as in Eq.(54). The resulting evolution is qualitatively similar to those in the 
Model 2 with Q„i < 0. Approximately, Pm{t) oc a(t)^^, Pr{t) oc a{t)~^, and Py{t) also has a scaling solution 
and exits the scaling regime, just like in the non-coupling case. The EoS of the YMC Wy approaches to —1 
from above, but does not cross over —1. To keep the paper short and concise, we will not repeat these details 
and not give the corresponding graphs here anymore. 

7. Asymptotic Behavior and Stable Attractor 

It is interesting to investigate the asymptotic behavior of the dynamical evolution. First we study the 
Model 1 with the YMC coupling to the matter component only. Now since the evolution of the the radiation 
component is independent of the the YMC and the matter, and the value of r at late time is much less than 
the other variables, so it can be neglected in the analysis of the fixed point. To find the fixed points, one sets 
dx/dN = dy/dN = in Eqs. (32) and (33), and obtains the relations at the fixed point: 

4 

Xf = —^yfe^^, (67) 



where the sub-index / refers to the the respective values at the fixed point. Prom these two equations one 
can write the asymptotic value of the fractional matter density 

^mf = (69) 

y/ - -3 

Eqs. (67) and (68) depend on the value of the ratio T/H, so does the solution (xf, yf). Here we consider 
the case that T/H is constant. Because i^mf must be larger than and smaller than 1, so yf must be in the 
range from —1 to at the fixed point. Thus, as long as F/H e (0, oo) there will exist fixed points. To be 
specific, consider T/H = 0.5. In the region of physics there is only one fixed point: 

(xf, yf) = (0.13666, -0.11499). (70) 

This is, in terms of the respective densities, 

(Pmf, Pyf) = ^6«;'(0.13666, 0.78888). (71) 

The stability of this fixed point can be analyzed in the conventional way as follows. Because the two equations 
for X and y are nonlinear, a local analysis can be given by linearizing the two evolution Eqs. (36) and (37). 
By a standard procedure, expanding x = Xf + e and y = yf + r], where e and rj are small perturbations 
around the fixed point, and keeping up to the first order of small perturbations, Eqs. (36) and (37) reduce 
to 

d I e 
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where M is a 2 x 2 matrix depending on the values of Xf, yf, and T/H, whose the elements are 
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where hf 



{2 + yfYhf_ 

(1 + yf)eyf + Xf + Vf. The general solution for the linear perturbations is of the form 

e = Cie^^^ + C2e^^^ (73) 

r, = Cae^^^ + de''^^ (74) 

where /xi and /X2 are eigenvalues of M. If they are both negative, the fixed point {xf,yf) is stable, and the 
solution is called an attractor. For the case of F/iJ = 0.5 one finds the matrix 



M 



-3.22146 
0.13180 



0.50112 
-2.17625 



and its two eigenvalues ni = —3.28123 and /X2 = —2.11648, respectively, both negative. Thus the fixed point 
of this model is stable, and is an attractor. For illustration, we plot in Fig. 16 the phase graph of trajectories, 
each trajectory starts with a different initial condition, and ends up at the fixed point of Eq.(70). As a matter 
of fact, one can check that the asymptotic behavior of other cases of Model 1 are also stable fixed points. 

The analysis of the asymptotic behavior can be done analogously for Model 3 with the couplings to the 
matter and the radiation. Setting dx/dN = dy/dN = dr/dN = in Eqs. (60), (61) and (62) one has the 
following three relations at the fixed point: 
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The fixed point {xf,yf,rf) as the solution of this set of equations depends on the ratios of rates T/H and 
T'/H as well. Consider the constant F and F'. The local analysis of the stability of the fixed point can 
be carried out similarly. By setting x = Xf + e, y = yf + tj, and r = r/ + 7, where e, j], and 7 are small 
perturbations around the fixed point, one has the equations 



d 

dN 







1 ^ 






1 


'1 




[1 



(78) 



where M' is a 3 x 3 matrix depending on the values ofxf,yf,jf,T/H and T'/H, whose the elements are: 
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Consider the specific case T/H = 0.5 and T'/H = 1.8 x 10"^. Substituting these into Eqs.(75), (76) and (77) 
yields the unique fixed point given by 



The three eigenvalues of the matrix M' are found to be —4.00015, —3.28124 and —2.11649, each being 
negative. Therefore, this attractor of the Model 3 is also stable. Notice that there are three quantities 
(.T, y, r) in the Model 3, so wc need the two phase graphs of trajectories. They are plotted in Fig. 17 for (x, y) 
and Fig. 18 for {r,y), separately. Each trajectory starts with a different initial condition, and ends up at the 
fixed point of Eq.(79). One can check the dynamics of other cases in Model 3 also have a stable attractor. 
As for Model 2 and Model 4, only Py{t) has an asymptotic constant value and has a stable attractor. 



Our motivation of this study is to investigate the coincidence problem for the cosmic dark energy in a 
spatially flat universe. We have presented a detailed and comprehensive analysis of the model of the effective 
YMC dark energy interacting with the matter and radiation. This work has been an extended development 
of our previous work on the non-coupling YMC dark energy model. Through the Priedmann equation and the 
dynamic equations for each cosmic component, once the couplings between these components are specified, 
the overall cosmic evolution is fully determined by the initial conditions of these three components. We 
have studied the evolution for the matter-dominated era starting from the equality of radiation-matter at 

7- ^ ^A^d 



{xf, yf, rf) = (0.13666, -0.11503, 4 x 10"^), 
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and the matrix 



/ -3.22149 0.50110 
0.13187 -2.17631 
^ -0.00008 0.00045 




8. Conclusion and Discussion. 



The major results of this work are the following. 

Given the initial dominant matter and radiation fi^j = flyi 2± 1/2 and the subdominant YMC energy 
density Clyi < 10"^, no matter what kind of coupling between the YMC and the matter, or between the 
YMC and the radiation, the evolution is such that the YMC is subdominant to, and keeps track of, the 
matter, until later, at a redshift z ~ 0.48 for a coupling > 0, or 2; ~ 0.35 for a coupling < 0, the 
YMC becomes dominant over the matter. The era is followed by a subsequent accelerating era driven by the 
dominant YM dark energy. As the evolution outcome, the universe arrives at the present state with ~ 0.7, 
Clr ~ 0.3, and Jl^ ^ W"^. It is very important to note that this has been achieved for a variety of coupling 
forms Qm and Qr, and, nevertheless, under a very broad range of initial condition Qyi ~ (IQ-^o, 3 x 10"^). 

If the YMC decays only into the matter, as a result of the coupling in Model 1, to achieve the present 
state of the universe, the decay rate needs to be of the same order of magnitude of the expansion rate of 
the universe, i.e., T ~ 0.5H. Moreover, for this coupled system, as t ^ 00, both py and pm asymptotically 
approach to constants, respectively. That is, for the system there is a unique attractor. Furthermore, as our 
analysis has shown, this attractor is stable. As an interesting behavior, the EoS for the YMC Wy always 
crosses over —1 around z ~ 2.5, and the present value is Wy ~ —1.1. This crossing —1 seems to be favored 
by the recent preliminary observations on SN la. 

When the YMC decays into both the matter and the radiation as in Model 3, there are two parameters 
r ~ 0.5H and F' ~ 1.8 x 10~'^H, representing the respective decay rate. The evolutional behavior is almost 
the same as Model 1, and Wy crosses over —1. Moreover, the radiation energy density also asymptotically 
approaches to a constant, ast^oo, like the YMC and the matter components. Most of the conclusions are 
the same as of Model 1. 

On the other hand, if the matter decays into the YMC as in Model 2 with a rate F ~ 0.02i7, or if both 
the matter and radiation decay into the YMC as in Model 4 with rates F ~ 0.0211 and F' ~ 1.8 x 10"^ H, 
respectively, then only py asymptotically approaches to a constant. Wy approaches to —1, but does not cross 
over —1. The evolution is nearly similar to the non-coupling case. 

Therefore, for all four types of models that we have studied, the coincidence problem can be naturally 
solved by introducing the effective YMC as the dark energy at the fixed parameter k given in Eq.(38). The 
present state of the universe is a natural result of the dynamic evolution. The past history of the evolving 
universe is that of the standard Big Bang model, and the future of the universe depends on the details of the 
coupling. If there is no coupling, or if the matter decays, or both matter and radiation decay, into the YMC, 
as in Model 2 and in Model 4, the matter and the radiation will keep on decreasing as Pmip) a{t)^'^ and 
pr{t) oc a{t)^'^. If the YMC decays into the matter only as in Model 1, then Pm{t) will asymptotically remain 
as constant, like Py{t) does, but Pr{t) oc a(t)~^. If the YMC decays into both the matter and radiation as in 
Model 3, then all the components Py{t), Pm{t), and Pr{t) will asymptotically remain as constant. In Model 
1 and Model 3 the future of the universe is a steady state, quite similar to that of the Steady State model, 
thus, in a sense, these two models bridge between the Big Bang model and the Steady State model. 

The distinguished characteristics of the YMC dark energy model are the following. 

The YM field is known to be indispensable to particle physics, the effective YMC employed in our work 
comes from quantum corrections up to 1-loop . Therefore, there is no room to adjust the form of the 



effective Lagrangian. This is in contrast to scalar field models, which have to design the form of potential 
and sometimes even the form of kinetic energy. 

The solution of coincidence problem has been relying on the parameter k in all our models. Viewed 
from the standard model of particle physics, the energy scale by k is much smaller than the other known 
microscopic energy scales. And this stands as the fine-tuning problem for any current cosmological model so 
far, and for our model as well. However, if the YM field in our model is regarded as a fundamental gauge 
field with k being the energy scale for this new physics, then the fine-tuning problem is traced up to the new 
physics. When the couplings are included, there are two more parameters T and V. But the present state of 
the universe requires that F be roughly the same order of magnitude of the expansion rate F ~ 0.5H, and F' 
be roughly three order lower. Since H ~ VGk? = K/nipi with lUpi being the Planck mass, thus the couplings 
F ~ K/nipi and F' ~ lO^'^K/nipi are also associated with the scale k. 

On the dynamic evolution, in comparison with scalar models, our models have the following features. All 
our models, for a broad range of the initial condition and for a variety of the coupling forms, automatically 
have the scaling property, i.e., Py{t) is initially subdominant to, and keeps track of to the matter. The 
accelerating stage begins only quite recently around a redshift z ~ (0.35, 0.48). This will allow the Big Bang 
cosmology to remain without drastic modifications. Besides, all our models have only one stable fixed point, 
uniquely determined by the ratio T/H and has nothing to do with the initial conditions ^lyi. Moreover, all 
the quantities in our model, especially Py{t) and Py{t) are continuous functions of t. So there is no Big Rip 
singularities in our models. Interestingly, as a function of t, Wy behaves quite smoothly during the evolution, 
going from ~ 1/3, approaching to —1. And in the case of the YMC decaying, Wy crosses over — 1 at z ~ 2, 
acquires the present value ~ —1-1, and settles down to an asymptotic value ~ —1.17. 
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Figure 1: Non-coupling case: The evolution of energy densities. For a wide range of initial conditions 
Pyi = (10~^°, 10~^)pmi; there always exists a scahng solution during the early stages, and Py{t) levels 
off and becomes dominant around z ~ 0.35. 

Figure 2: Non-coupling case: The evolution of fraction of energy densities. The two curves of Q,y for 
fji — l and iji — 16.12 are very close to each other, and nearly overlap in this figure. 

Figure 3: Non-couphng case: The evolution of EoS. Wy approaches to — 1 as i ^ oo, but does not 
cross over —1. 

Figure 4: Model 1 with Qm > for T/H = 0.5 : The evolution of energy densities with the YMC 
decaying into the matter. For a wide range of initial conditions Pyi = (10^^", 10^'^)pmi, i-e., yi = 
(1, 15), there always exists a scaling solution during the early stages, and Py{t) levels off and becomes 
dominant around z ~ 0.48. Due to the coupling, Pm{t) also levels off at late time. 

Figure 5: Model 1 with Qm > ioi T/H — 0.5 : The evolution of fractional energy densities in the 
same model as in Fig.4. 

Figure 6: Model 1 with > ior T / H — 0.5 : The evolution of EoS in the same model as in Fig.4 
and Fig. 5. Due to the coupling, Wy crosses over —1, takes on value ~ —1.1 at 2; = 0, and approaches 
—1.17 asymptotically. This is in contrast to the non-coupling case. 

Figure 7: Model 1 with Qm > for F/H ^ 0.25ef f±| : The evolution of energy densities with the 
YMC decaying into the matter. For a wide range of initial conditions pyi = (10~^°, 10~^)pmj, there 
always exists a scaling solution during the early stages, and Py{t) levels off and becomes dominant 
around z ~ 0.48. Due to the coupling, pm{t) also levels off at late time. This is quite similar to the 
case of a constant rate T/H = 0.5 in Fig.4. 

Figure 8: . Model 1 with Qm > for T/H = 0.25e^ fj| : The evolution of fractional energy densities 
in the same model as in Fig. 7. 

Figure 9: Model 1 with > for F/iJ = 0.25e^ 1^ '■ The evolution of EoS in the same model as in 
Fig. 7 and Fig. 8. Again, due to the coupling, Wy crosses over —1, and takes on value ~ —1.1 aX z = 0. 

Figure 10: Model 2 with < for the constant coupling T/H — 0.02 : The evolution of energy 
densities with the matter decaying into the YMC. 

Figure 11: Model 2 with Qm < for the constant coupling F/i? = 0.02 : The evolution of fractional 
energy densities in the same model as in Fig. 10. 

Figure 12: Model 2 with Qm < for the constant coupling T/H = 0.02 : The evolution of EoS in the 
same model as in Fig. 10 and Fig.ll. Since the YMC gets energy from the matter, Wy approaches to 
— 1, but does not cross over —1. 



Figure 13: Model 3 with Qm > and Qr > for F/H = 0.5 and F'/H = 0.00018 : The evolution of 

energy densities with the YMC decaying into both the matter and the radiation. For a wide range of 
initial conditions Pyi = (10"^'^, 10^^) Pmi, there always exists a scaling solution during the early stages, 
and Py{t) levels off and becomes dominant around z ~ 0.48. Note that, due to coupling, both Pm{t) 
and pr{t) level off like Py{t). 



Figure 14: The evolution of fraction of energy densities in the same model as in Fig. 13. 



Figure 15: Model 3 with Qm > and Qr > for F/H = 0.5 and F'/H = 0.00018 : The evolution 
of EoS in the same model as in Fig. 13 and Fig.l4. Note that Wy also crosses over —1 and takes on a 
value ~ —1.1 at 2; = 0. 



Figure 16: Model 1 with Q„i > for F/H = 0.5 : The trajectories in the phase plane. Each 
trajectory starts with a different initial condition. All of them approach the fixed point {xf,yf) — 
(0.13666, -0.11499). 



Figure 17: Model 3 with Q„ > and Q, > for F/i/ = 0.5 and F'/H = 0.00018 : The trajectories in 
the phase plane {x,y). Each trajectory starts with a different initial condition. All of them approach 
the fixed point {xf, yf, r/) = (0.13666, — 0.11503, 4 x 10~^). The parameters are the same as in 
Fig.l3. 



Figure 18: Model 3 with Qm > and > for F/H = 0.5 and F'/H = 0.00018 : The trajectories 
in the phase plane {r,y). The parameters are the same as in Fig. 13. 
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